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Abstract 

Using 5.2 fb~^ of e~^e~ annihilation data recorded with the BABAR detector at the PEP-II storage 
ring while operating on the T{4S) resonance, a sample of fully reconstructed decays in the 
hadronic modes D*~tt~^ and B^ D*~p~^ have been reconstructed. In this paper, a 

study of these events is reported, including preliminary measurements of the absolute branching 
fractions for these modes, which are found to be B{B'^ D*^tt^) = (2.9 ± 0.3 ± 0.3) x 10~^ and 
e(B° D*-p+) = (11.2 ± 1.1 ± 2.5) X 10-3. 
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1 Introduction 



The dominant hadronic decay modes of the B meson leading to open charm in the final state 
involve tree-level diagrams where the b ^ c transition leads to a charmed meson and an external 
W, which often emerges as a charged vr, p, or ai. In this study, we report a new measurement of 
the branching fractions for D*^tt^ and B^ D*^ Q The former is known to about 10% 

from previous measurements, while the latter has proven more difficult to determine due to the 
helicity of the p |, |, |]. 

2 The BaBAR detector and dataset 

The data used in this analysis were collected with the BABAR detector while operating in the PEP-II 
storage ring at the Stanford Linear Accelerator Center. For this analysis, we use a sample equivalent 
to 5.22 ±0.16 fb~^ of integrated luminosity, collected while running on the T(4S') resonance. This 
corresponds to (5.93 ± 0.21) x 10^ BB pairs, assuming the T(45) decays only to B mesons. 

The BABAR detector is described in detail elsewhere [Q]. Charged particle tracking is provided 
by a 5-layer silicon microstrip detector, capable of stand-alone pattern recognition, and a 40-layer 
cylindrical drift chamber. This is followed outward in radius by a specialized particle identification 
system, based on detection of Cherenkov light generated in quartz. A calorimeter consisting of 
thallium-doped cesium iodide crystals provides detection for photons, and particle identification 
for electrons. These devices sit inside the superconducting coil, which provides a 1.5 T field with 
3% uniformity over the tracking volume. The flux return for the magnet is finely segmented and 
instrumented with resistive-plate chambers to provide both muon identification and crude hadronic 
calorimetry for the detection of mesons. 

3 B reconstruction method 

Hadronic events for this study are selected by a set of simple requirements designed to minimize 
systematic errors in B counting. We require more than three charged tracks forming a primary 
vertex within 0.5 cm of the beam spot in both transverse directions, a sum of charged and neutral 
energy greater than 5 GeV, and the normalized second Fox- Wolfram moment for the event, R2, 
calculated from charged tracks and neutrals in the T(45) frame, less than 0.5. Neutrals are clusters 
in the calorimeter with no associated tracks, more than 30MeV in energy, and a lateral moment 
of the shower distribution less than 0.8. The R2 requirement is designed to reject the jet-like 
continuum events over the more uniformly distributed T(45) decays. A fiducial requirement of 
0.410 < 9Lab < 2.540 for charged tracks and 0.410 < Oiab < 2.409 for neutrals, where Oiab is the 
polar angle to the beam line, is made in calculating these quantities and in the subsequent steps 
of the study. 

B^ candidates in the channels D*^tt^ and D*^p~^ are reconstructed using the mode D*~ — > 
D^TT^, followed by —>■ K~^7r~. The is seen in its decay to vr'^'Tr'^. All charged tracks are 
required to originate close to the beam spot; the distance of closest approach in the transverse 
plane, is required to be less than 1.5 cm and along the beam line to be less than 3 cm. The 
daughters of the must have a transverse momentum, pT, greater than 100 MeV/c, and include 

^Here, and throughout this document, we use the convention that a particular candidate state also implies the 
charge conjugate state is included. 
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at least 20 drift chamber hits. The soft pion track is only required to havepx greater than 70 MeV/c, 
taking advantage of stand-alone track finding in the silicon detector to improve acceptance. No 
particle identification is used for this analysis. Neutral pions are formed from pairs of photons with 
energy greater than 30MeV and having a lateral shower moment less than 0.8. The invariant mass 
of the candidate must lie within ib25MeV/c'^ of the nominal vr'^ mass. 

candidates are formed from K~^7r~ combinations, where the kaon and pion are required to 
have a minimum momentum of 200 MeV/ c, and the invariant mass of the candidate must be within 
±2.5(T of the nominal mass [^. All candidates are required to have momentum greater than 
1.3GeV/c in the T(4S') frame. A vertex constraint fit is performed, for which the probability 
must be greater than 0.1%. 

The reconstructed is then combined with a soft pion having charge opposite to the kaon 
to form D*~ candidates. A vertex constraint fit is applied to improve the angular resolution for 
the soft pion, using a fixed effective vertical beam spot size of 40 fim. In those cases where the 
fit converges, D*~^ candidates are selected by the requirement that Am = m{D^-K~) — m{D^) lies 
within ±2.5(7 of the nominal mass difference The width is taken from a weighted average of 
the two-Gaussian distribution that is required to fit the Am distribution. 

_gO _^ D*-TT+ candidates are reconstructed by combining a D* and a vr^ with momentum 
greater than 500MeV/c. The decay D*~7r~^ involves a pseudoscalar initial-state particle 

decaying into a vector and pseudoscalar, so that the final-state D*^ is polarized. Therefore, the 
helicity angle 9h between the soft pion direction and the D*~ direction in the D*~ rest frame should 
be distributed as cos^ Oh- In contrast, combinatorial background is uniform in cos 6h- Therefore, 
B meson candidates are selected with the additional requirement | cos^^^l > 0.4. This requirement 
removes 40% of the background and only 5% of signal events. 

For the B^ — > D*~ mode, p'^ candidates are formed by combining a vr'' meson and a charged 
pion with momentum greater than 200MeV/c. We require the momentum of the p"*" to be greater 
than IGeV/c and the tt^tt^ mass to lie within 150 MeV/ of the Z?"*" mass. Finally, the opening 
angle, 9t, between the thrust angle of the B candidate and the thrust axis of the remaining 
charged and neutrals in the event is required to satisfy | cos 6t\ < 0.9 in order to remove continuum 
backgrounds. 

In the case of a correctly reconstructed B meson produced by the decay of an Y{4S), within 
the experimental resolution, the measured sum of neutral and charged energies, E^, must be equal 
to the beam energy, E^, both evaluated in the T{4S) frame. We define AE to be the difference 
between the measured B candidate energy and beam energy in the T(4S') frame, AE = E^ — E^. 
The resolution on ctab is predicted to be 16.0 it 0.6 MeV for D*~tt~^ and 28 MeV for D*~ . The 
beam-energy substituted mass, m-Es, is defined as m'^g = (E^)'^ — {J2iPi)'^i where the Pj is the 
momentum of the ith daughter of the B candidate. The predicted resolution in mss is typically 
about 2.5 MeV/c^ for D*-7r+ and 3.1 MeV/c^ for D*' p+ . This is about a factor of 10 better than 
the resolution in the reconstructed invariant mass. The resolution for mss is dominated by the 
beam energy spread rather than by the detector resolution. It is largely uncorrelated with the error 
on AE. 

The variables AE and tties are used to define a signal region and also sideband regions for 
background studies. For all modes, the region between 5.2 and 5.3 GeV/c^ in rriES and between 
±300 MeV in AE is used to study the B candidates. The peak position, m^, which should be 
the nominal B mass, and the resolution (Jm-p^Q, are extracted from the distribution of mss after 
requiring AE be consistent with zero to within ±2.5(t. The resolution in AE is extracted from the 
AE distribution obtained by requiring that mss lie within ±2.5o"r7iEs "^B- 
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The signal region in the two dimensional plane m-^s versus AE is defined ±2.5(7 wide 

centered at the nominal B mass, m^, and at AE = 0. By staying below [A^'l = ttt-tt, we avoid 
correlated background from B decays where a real final-state charged pion is either not included 
in the reconstruction or a random soft pion from the recoil B is added to the observed state. 

Only one candidate per event is allowed to appear in the mES versus AE distribution. The 
criteria selected is to consider only the entry with the smallest absolute value for AE. The resulting 
two dimensional distribution of candidate events in AE and rriEs is shown in Fig. |^. 




mp„ (GeV/c") 




mp„ (GeV/c") 



Figure 1: Distribution in data of AE versus tues for B candidates in the channels (a) B^ — > D* ir^ 
and (b) S° ^ D*-p+. 



4 Branching fraction measurements 

The measurement of branching fractions requires an estimate of the combinatorial background in 
the TTT-Es distribution near the nominal B mass. For the channel B^ D*~tt~^ , the m-Es distribution, 
shown in Fig. is fitted with a background function Q which parameterizes how the phase space 
approaches zero as the energy approaches E^: 



fBG{i) 



Nidi-e^M^'^-n) 



(1) 



where ^ = niEs/El, and the normalization and the shape are determined by the parameters N and 
K. The signal is characterized by a Gaussian distribution with free mass, m^, and width, Crn-Es- "^^^ 
projection of the signal as a function of AE, also shown in Fig. is fitted using a linear background 
function plus a single Gaussian distribution with free mean and width, cr a_b- The region —0.3 < 
AE < — O.lSGeV, containing feeddown from D*~ p where a pion is left out of the reconstruction, 
is excluded from the fit. The observed width of the signal in mEs is 2.45 ± 0.18MeV/c2 while the 
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Figure 2: Distribution of m^s ^or \AE \ < 2.5<ta_b (left), and AE for jm-Es — < 2.5am-Es 
(right) in the channel — > D*~'ir^. The fits are described in the text. 

Table 1: Observed and expected yields and efficiencies for decay modes. 



B^ mode 


Observed yield 


Efficiency (%) 


D*-TT+ 

D*-p+ 


119 ± 11 
131 ± 13 


27.0 ± 1.0 
7.6 ±0.6 



AE resolution is observed to be 18.8 it 0.9 MeV. Acceptances are calculated in terms of the fitted 
widths observed in data and Monte Carlo simulation. 

For the channel, B^ D*~ p'^ the observed m^s distribution is shown in Fig. ^. This has been 
fitted with the same background function and signal function described above. The projection of 
the signal as a function of AE, also shown in Fig. ^, is fitted using a linear background function 
plus a single Gaussian distribution with free mean and width, (Tab, with the region —0.3 < AE < 
—0.13 GeV once more excluded. The width of the signal in rriEs is 3.5 ± 0.3 MeV/ while the AE 
resolution is observed to be 39.5 ± 4.7 MeV. 

A detailed Monte Carlo simulation is used to determine the acceptance for the signal events. 
Control samples are used to determine the uncertainties on crucial performance characteristics of 
the simulation. The fitted yields for signal events and estimated efficiencies for each channel are 
listed in Table |^. 

A variety of sources contribute to the systematic errors on the final branching fraction results. 
The number of produced B mesons is extracted from the ratio of multihadron-to-muon pair events 
on- and off-resonance, after application of the simple event selection criteria described above. Af- 
ter extrapolation from accepted to produced numbers of BB pairs, the estimated uncertainty is 
3.6% §. 
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Figure 3: Distribution of m^s ^or \AE \ < 2.5(7^e (left), and AE for |mEs — < 2.5am-Es 
(right) in the channel — > D*~p^. The fits are described in the text. 

The primary check on charged track efficiency is obtained by studying the probabihty for ob- 
serving drift chamber versus sihcon detector-only tracks in inclusive — > i^^vr+vr+Tr". This is 
compared with the rate for finding the third track in one-versus-three topology tau-pair decays. A 
final check is the observed multiplicity distribution in T{4S) events. A systematic error of 2.5% 
per track with pT > 1 GeV/ c on the overall efficiency scale is determined by comparing the three 
methods. The soft pion efficiency is determined from a study of the forward-backward asymmetry 
in observed D*~ decays, and is also constrained to some extent by the same charm and tau-pair 
studies. The pT resolution in the Monte Carlo events is adjusted to reproduce the resolution seen in 
cosmic ray muons. The overall systematic uncertainty in the final result due to tracking efficiency 
is estimated to be 7.9% for the final states considered in this study. An additional systematic error 
contribution come from changes in the ratio of observed and predicted efficiencies as the selection 
requirements are varied within a reasonable range (3%). 

In the case of the — > D*~ channel, there are additional systematic errors due to the 
p"*" and its polarization. The reconstruction efficiency for the vr*^, as modeled in our Monte Carlo 
simulation, is verified to be accurate to within 5% Q. The width of the signal in the AE projection 
is modestly larger than predicted by the Monte Carlo simulation. Therefore, we have varied the 
AE requirement for the signal band over a wide range around the nominal 2.5it, and estimate a 
further 10% uncertainty due to this requirement. Finally, we have extracted the branching fraction 
separately in the two hemispheres of the p helicity angle, 9 nip), greater and less than zero. A further 
15% systematic error is assigned on the basis of the observed difference in branching fraction. The 
helicity distributions for the p and the D*~ are consistent with expectations based on previous 
measurements and our Monte Carlo efficiency calculation, although we have not yet attempted 
to determine the polarization from the data. 

An additional check of the result for B^ — > D*~7:~^ is to relax the cos 9h requirement, and 
extract the fitted signal as a function of cos9h instead. The distribution of the B candidates 
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Figure 4: Distribution of D* vr^ as a function of the helicity angle of the soft pion, cos 9h, 

in the D*~ rest frame after acceptance correction. The vertical axis has arbitrary units. 

within ib2.5fT772gg of the nominal B mass and ±2.5cr/\E of AE = exhibits no forward-backward 
asymmetry after acceptance corrections as can be seen in Figure This confirms our understanding 
of the soft pion efficiency in the D*~ decay. Overlayed on the 126 signal candidates, with an 
estimated background of 4 events, is a fit with cos'^ 9h, yielding a for goodness-of-fit of 8.6 for 
7 degrees of freedom. 

The PDG compilation |^ of measured branching fractions for the D*^ D'^tt'^ , (68.3 it 1.4)%, 
and K~Tr~^, (3.83 it 0.09)%, are used in computing our final results. The measurement 

errors on these branching fractions are included in the systematic error on our final result. We also 
assume the T(4S') decays into B^B^ pairs with a 50% fraction; no systematic error is assigned to 
this value. 

Based on fitted yield of signal events, the estimated efficiency, and the number of produced B 
mesons in our sample, the preliminary results for the branching fractions for B^ D*~7r^ and 
^0 _^ D*~p+ are (2.9 ± 0.3 ± 0.3) x 10"^ and (11.2 ± 1.1 ± 2.5) x lO"^ respectively. The branching 
fraction for B^ — > D*~ includes all non-resonant and quasi- two-body contributions that lead to 
a ir'^ir^ invariant mass in the p band. However, the acceptance for non-resonant D*~7:~^7:^ decays is 
about 15% of D*^ p^ so that, combined with the known branching fraction for this mode, the non- 
resonant contribution to our result for B^ D*^ p'^ is expected to be quite small. Both branching 
fraction results compare well with previous measurements and with the world average |^. 

5 Summary 

B^ decays to D*^tt^ and D*~ p'^ have been studied using the decay chain D*~ — > D^'it~ , followed 
by — > K~^7:~ . The preliminary branching fractions obtained for these channels, (2.9ib0.3ib0.3) x 
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10 and (11.2 it 1.1 it 2.5) x 10 respectively, are compatible with previous observations [0, 
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